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Abstract 

The isotropic XY model (s = 1/2) in a transverse field, with uniform long-range interactions 
among the transverse components of the spins, on the inhomogeneous periodic chain, is studied. 
The model, composed of N segments with n different exchange interactions and magnetic moments, 
is exactly solved by introducing the integral gaussian transformation and the generalized Jordan- 
Wigner transformation, which reduce the problem to the diagonalization of a finite matrix of nth 
order. The quantum transitions induced by the transverse field are determined by analyzing the 
induced magnetization of the cell and the equation of state. The phase diagrams for the quantum 
transitions, in the space generated by the transverse field and the interaction parameters, are 
presented. As expected, the model presents multiple, first- and second-order quantum transitions 
induced by the transverse field, and it corresponds to an extension of the models recently considered 
by the authors. Detailed results are also presented, at T = 0, for the induced magnetization and 
isothermal susceptibility Xt ^ function of the transverse field. 
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I. INTRODUCTION 



The study of the critical quantum behavior of systems[l|], which is induced by quantum 
fluctuations, has been object of great interest in recent years. This critical behavior, which 
controls the properties of the systems at very low temperature, is present in different sys- 
tems. In particular, for magnetic systems, they have been responsible for unusual properties 
observed in low dimensional magnetic materials a a. Therefore, the study of the criti- 
cal behavior of spin systems in low dimension, particularly the exactly soluble, is of great 
importance for understanding the properties of these materials. 

Among these models, the one-dimensional XY model introduced by Lieb, Schultz and 
Mattis [4] , despite being almost fifty year old, is still the best one to describe exactly magnetic 
quantum transitions. The rather rich quantum critical behavior presented by the model can 
be seen in a recent work by the authors [sl and in the references therein, where they study 
the anisotropic model on the inhomogeneous periodic lattice. The study of the model on 
the inhomogeneous open lattice has also been recently addressed by Feldman 61]. 

The isotropic model on the inhomogeneous lattice has also been studied by the authors [3] , 
where a detailed study of the static and dynamic critical properties is presented. Although 
the model has been applied mostly in the study of the quantum critical behavior of magnetic 



systems 



0, y, y, y, 



10|, more recently it has also been applied in the study of quantum 



entanglement, which plays an essential role in the quantum computation. These applications 
can be found in the recent work by Amico et a/. and in the references therein. In 
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particular, in a recent work on quantum communication in a spin system, Avellino et al. 
have studied the strong effect of long-range interaction on the fidelity of transmission of 
quantum information. 

Besides the importance of the long-range interaction in the transmission of quantum 
information in spin chains, its presence can induce classical critical behavior in these systems 
which is essential for the study of the classical/quantum crossover. 

The one- dimensional XY model is among the models which present this behavior, pro- 
vided a uniform long-range interaction is considered. In particular, for the isotropic one, in 
the presence of a homogeneous long-range interaction along the transverse field direction, 
it can still be solved exactly, and its solution has been obtained by the authors [l^ . 



Besides the appearance of classical critical behavior, the most important features presented 
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by the homogeneous model have been the existence of quantum bicritical points and, par- 
ticularly, the existence of first order quantum transitions. Another relevant result obtained 
in the study of the homogeneous case was the exact determination of the classical/quantum 
crossover for first and second order phase transitions, and the verification of the scaling 
relations proposed Continentino and Ferreira |il5i] for first order quantum transitions. 

As pointed out by Pfieiderer 16|, quantum first order transitions can be driven by dif- 
ferent mechanisms, and to look at this new quantum critical behavior, in an exactly soluble 
problem, has been the one of the motivations to analyze the isotropic model on the inho- 
mogeneous periodic chain with long-range interaction. Besides this point, the possibility of 
having multiple first and second order quantum transitions and multicritical points of higher 
order have also been important motivations for considering the long-range interaction in the 
inhomogeneous model. 

Therefore, in this paper we will extend our previous results [3, Q, [l^ by looking at an 
extended version of the models previously considered. Although the model presents classical 
and quantum critical behavior, we will restrict our analysis to the quantum one. 

In Section 2 we present the exact solution of the model and obtain the functional of the 
Helmholtz free energy at arbitrary temperatures. The quantum critical behavior is discussed 
in Section 3, and it is determined from the equation of state at T = 0. Explicit results for 
the quantum phase diagrams, induced magnetization and isothermal susceptibility Xt 
for systems with different sizes of the unit cell are also presented. Finally in Section 4 we 
summarize the main results of the paper. 



II. THE MODEL AND THE FUNCTIONAL OF THE FREE ENERGY 

We consider the isotropic XY model (s = 1/2) on the inhomogeneous periodic chain 
with N cells, n sites per cell, and lattice parameter a, in a transverse field, with long- 
range interactions among the spin components in the z direction. The unit cell of the 
inhomogeneous lattice is shown in Fig. 1, and the Hamiltonian of the model in its more 
general form is given by 



4 



d=na 



• m 







m=l m=2 



m=n 



l-th cell 



FIG. 1: Unit cell of the inhomogeneous chain. 
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where the parameters Jm are the exchange coupling between nearest-neighbours, fim the 
magnetic moments, h the external field, J' the uniform long-range interaction among z 
components, and where we have assumed periodic boundary conditions. If we introduce the 
ladder operators 

S^ = S'=± iS\ (2) 



and the Jordan- Wigner transformations 



»n— 1 



l'=l [^m'=l 

-r^ r, -- 
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tonian as 



and C; ^ are fermion annihilation and creation operators, we can write the Hamil- 
17] 
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(6) 
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and 



I±P 



(7) 



N n 



with P given by 

P = exp imJ2Yl ^lmCi,m 1 ■ (8) 

As it is well known |17l . Il8l . Il9l |. since the operator P commutes with the Hamiltonian, 

the eigenstates have definite parity, and P~(P+) corresponds to a projector into a state of 
odd (even) parity. 

Introducing periodic and anti-periodic boundary conditions on c's for and H^, re- 
spectively, the wave-vectors in the Fourier transforms 



(9) 



and 



1 ^ 

Ag^^ = -j= ^ exp {iqdl) 
^ ^ 1=1 



_are given by g =77^ for periodic condition and = ^^^77^, for anti-periodic condition 



20|], with I = 0, ±1, ±N/2, and H and can be written in the form 



(10) 



where 

n 



m=l 
n-1 



m=l 
J, 



J A J± ^^v4g±,i exp (-2cig^) + -^v4,±,„ exp(zdg^) 



+ 



+ 7, ,^Aj±^,exp(-^c/g±) + A,±,i V,n exp(«c/g±))] - ^ (^J±,™V 



m=l 



(11) 



Although H and if"'" do not commute, it can be shown that in the thermodynamic limit 
all the static properties of the system can be obtained in terms of H~ or |17l . llSl . Il9 |. 
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Therefore, by considering periodic boundary conditions on c's, we can identify H = H and 
q = q" , Hq can be written in the form 

m=l '- \ / J r?i=l 



- y [^J,n^?,iexp(-i(ig) + Aj iAg,nexp(idg) 
The partition function is then given by 
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J' 



q,m 



\m=l 



■IN 




n (^^) ' 



(12) 
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(14) 
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Al^Aq^i exp(-idg) + Aj,iAg,n exp(idg)] + (/x^/i - J') A\ 



m=l 



(15) 



m=l 



Since the long-range interaction term commutes with the Hamihonian, we can introduce the 
gaussian transformation 

oo 

2 



exp (62) 



' exp ( h \/26a; ] dx, 



(16) 



in the previous expression, so that the partition function can be rewritten in an integral 
representation as 



Zg = exp I ( ^'"^ ~ y 



m=l 



exp I -y 1 Tr 



exp -PHq{x) 



dx, 



where x = x/\/nN, and the effective Hamiltonian Hq[x) is given by 
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Hq{x) = 
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By introducing the canonical transformations 

n n 
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and by imposing the condition 



[^q,k, Hq{x)] — Sq,kiq,k, 



this leads, for the coefficients Uq^km, to the equation 
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and the u's satisfy the orthogonality relations, 
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(20) 
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(22) 
(23) 



and hj — — J' + y^x. Therefore the effective Hamiltonian can be written in the 
diagonal form 

Hq{x) = J2^1'''^l,kVq,k: (24) 
k 

where the spectrum Sg of Hq is determined from the determinantal equation 



det{Aq - EqT) = 0, 



(25) 
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and the operators ri's are given in terms of A's and c's by the expression 

N n 



Vq,k = K,km\m = i Yl '^^^^^^ <,fcmQ,m, (26) 



m=l * 1=1 m=l 



which have been obtained by using eqs.Q and fl22|) . From this result the partition function 
can be written in the form 



Z 



N 




exp 



( J' 



m=l 
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exp In (1 + exp{-Peg^k))^ dx. 



(27) 
(28) 



In the thermodynamic hmit, the partition function can be evaluated by Laplace 's method 
by imposing the condition that g'ixo) = 0, where g{x) is given by 



g,k 



and xq is equal to 



Xo 



-E- 

nN ^ 1 

q,k 



VWJ" 



+ exp(— /3e 
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which can be written in terms of the average induced magnetization M^, defined as 

= ~^Y^^"'^^t"^^ ^ ~^Y1 ^^rnUlkmUq,km{vlkVg,k) - ^, 
l,m q,k,m 



(29) 



in the form 



y^2pj 2 



Finally, we can obtain from the previous results the functional of the Helmholtz free 
energy per lattice site which is given by 



nN 



(30) 



q,k 



The equation of state is obtained numerically from this functional by imposing the con- 
ditions 



•' 0, with — ^ > 0. 



dM 



(31) 
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The numerical solution is more easily obtained for n < 4, where there are analytical solutions 
for eqs. fl20l) and fl25|) [7], and in the cases where we have uniform magnetic moments, namely, 
fijn = fJ'- In this situation, the term — X]"™ f^^^im commutes with the Hamiltonian, and 
consequently the effect of the field is to shift the spectrum. This means that the solution of 
eq. (|2^ can be written as 

1 2 V 

e,,k = -fih + J'-J —x-o = e^^, -fih- 2 J'M^ (32) 



where ?^ ^ is the energy of the excitations of the model for zero transverse field and 
interaction, and we identify the expression fih + 2J'M^ as an effective field hejf 



ong-range 
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III. THE QUANTUM CRITICAL BEHAVIOR 

In the limit T — 0, the functional of the Helmholtz free energy per lattice site, eq. (l30!l . 
can be explicitly written as 

f = ^ = ^ + J'M^(M^ + 1) - —V / e,,,dq,, (33) 
iV / nn ^ — ^ / 

where g^, is obtained by imposing the condition Sq^k = and we have considered the lattice 
spacing a = 1 and uniform magnetic moments /i = 1. 

Under these conditions the average induced magnetization, M^, given in eq. fl2^ . can also 
be written in the form 

= V [ sign{eg,k)dqk. (34) 
Znn^ — 'J 



k=l 







The quantum phase diagram is determined from the equation of state, which is obtained 
numerically from the previous equations by considering the conditions shown in eq. fl3ip . and 
in particular we can determine the magnetization as a function of the transverse field. As 
it is well known, the second order phase transitions are determined by imposing the limit 
— >■ and it can be shown numerically, for arbitrary unit cell sizes, that the critical field 
varies linearly with the long-range range interaction J' and that these phase transitions, as 



in the homogenous model IJ], occur for J' < only. 



The first order phase transitions are determined by imposing the additional condition 
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/(M^) = / (M;) , (35) 

where are the magnetization plateaus, which are identical to the ones in the model 
without the long-range interaction |7^. This is due to the fact that, for uniform magnetic 
moments, the long-range interaction does preserve the azimutal symmetry of the model and 
consequently the magnetization plateaus satisfy the quantization condition j2l| 



n 



I - ) = /i X integer, (36) 



which does only depend on the symmetry of the Hamiltonian. 

As in the homogenous case |[14j], the first order phase transitions occur for J' > only, 
and in this case, besides the bicritical points where the second order lines meet the first 



order ones ISj, there are triple points which correspond to the point where three first order 
lines meet. 

From the magnetization as a function of the field we can obtain the isothermal suscepti- 
bility, Xt^ "which is given by 

XT - (37) 

n oh 

The main results are shown in Figs. 2 to 11. In Fig. 2 we present the phase diagram for 
= 2, and as expected there are two second order lines for J' < which end up at bicritical 
points for J' = 0. For J' > 0, there are two first order lines which collapse into a single one 
at a triple point. In Fig. 3(a) we present the magnetization as function of the field, and as 
it can be seen the different critical behaviors are explicitly shown depending on the sign of 
the long-range interaction. It should be noted that as we approach the triple point the two 
first order phase transitions collapse into a single one, which corresponds in this case to a 
jump in the induced magnetization from zero to one half. 

As in the homogeneous model|l4|, there is a universal curve to which all magnetization 
data collapse, independently of the order of the transition. This is presented in Fig. 3(b), 
where we show the magnetization as a function of the effective field h^ff = h + 2J'M^. 

In Fig. 4 we present the quantum phase diagram for = 3. As in the previous case, 
for J' < we have second order phase transitions and for J' > we have first order 
phase transitions. There are three second order lines and two first order lines meet at a 
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unique triple point. The magnetization is presented in Fig.5(a) for different values of J', 
which characterize the different behaviors, and in Fig 5(b) we present the collapse of the 
magnetization when plotted as a function of the effective field /ie//- 

It should be noted that for a different set of parameters we could have two triple points, 
and in this case we will have a single first order line beyond the critical value J', which is 
associated to the second triple point. This situation can be seen in the phase diagram shown 
in Fig 6, for n = A, where it is explicitly shown that we can have three triple points. Since 
we have no suppression of a phase^ransition, there are four second order transition lines. 
As it has been shown for J' = [7], an adequate choice of the exchange parameters can 
suppress a phase transition. For n = 4, this condition corresponds to J1J3 = J2J4 and the 
phase diagram for this case is presented in Fig. 7. As it can be verified in this figure, we 
have three second order lines, instead of four, and just one triple point. 

The magnetization associated to the parameters defined in Fig. 6 is presented in Fig. 8. 
As in the previous cases, in Fig 8(a) we have the magnetization as a function of the field 
for different values of J', and in Fig. 8(b), the universal curve for the magnetization as a 
function of the effective field heff. In Fig. 9, we present the phase diagram for n = 5. Since 
n is odd, there is no suppression of any transition for J' = j3] , and, as expected, there are 
five second order transition lines. The associated magnetization as a function of the field is 
presented in Fig 10(a) and the magnetization universal curve is shown in Fig. 10(b). 

Finally, in Fig. 11, we present the isothermal susceptibility, Xt^ n = 3 and different 
values of J'. As in the homogeneous model, the isothermal susceptibility, at the second order 
phase transitions, diverges for J' = only, and its multiple phase transitions have the same 
critical exponents of the homogeneous model and consequently belong to same universality 
class. 

IV. CONCLUSIONS 

In this work we have considered the isotropic XY model with an uniform long-range 
interaction along the z direction, in a periodic inhomogenous lattice with N cells and n sites 
per cell. The exact solution of the model was formally obtained at arbitrary temperatures 
and distribution of magnetic moments and exchange constants. Explicit equations have 
been obtained for the functional of the Helmholtz free energy from which the equation of 



12 



state can be determined numerically. 

The analysis of the critical behavior has been restricted to the quantum phase transitions 
and we have shown, as in the homogeneous model, that the system presents first order 
phase transitions when the long-range interaction is ferromagnetic, and second order phase 
transitions when the long-range interaction is antiferromagnetic. The model also presents 
multiple first and second order phase transitions and the number of critical lines is equal to 
ra, for n odd, and less than n, for n even, provided the exchange constants in the xy plane 
satisfy a special relation. 

The second order critical lines meet the first order ones at J' = 0, which are bicritical 
points and there are multiple triple points, where three first order critical lines meet, de- 
pending on the unit cell size. The critical exponents have been obtained numerically, and it 
has been shown, as expected, that the model belongs to the same universality class of the 
homogeneous one. 

Finally, we would like to point out that is of paramount importance the presence of mul- 
tiple first order quantum transitions and triple points, which we have shown exactly to exist 
in the model, since we believe that this behavior is related to the different mechanisms from 
which the first order phase transitions are driven. These mecanisms have been thoroughly 
discusssed by Pfieiderer[16|, by analyzing quantum critical behavior obtained for different 
materials. 
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FIG. 2: Phase diagram for the quantum transitions as a function of the strength of the long-range 
interaction J' /Ji, for n = 2 and Ji = 1, J2 = 2. For J' / J\ > 0, the critical lines identify the 
first-order phase transitions and, for J'/Ji < 0, the second order phase transitions. 
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FIG. 3: (a) Magnetization as a function of /i/Ji, and (b) universal curve for the magnetization 
as a function of the effective field heff/Ji (heff = h + 2J'M^), at T = 0, for different values of 
r(r = J'/Ji), in the regions where the system undergoes first (r > 0) and second-order (r ^ 0) 
quantum transitions, for n = 2 and Ji = 1, J2 = 2. 



16 




-2-1012 
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FIG. 4: Phase diagram for the quantum transitions as a function of the strength of the long-range 
interaction J' /Ji, for n = 3 and Ji = 1, J2 = 1, J3 = 2. For J' / J\ > 0, the critical lines identify 
the first-order phase transitions and, for J'/Ji < 0, the second order phase transitions. 
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FIG. 5: (a) Magnetization as a function of h/Ji, and (b) universal curve for the magnetization 
as a function of the effective field hgff/Ji (heff = h + 2J'M^), at T = 0, for different values of 
r(r = J'/Ji), in the regions where the system undergoes first (r > 0) and second-order (r ^ 0) 
quantum transitions, for n = 3 and Ji = 1, J2 = 1, J3 = 2. 
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FIG. 6: Phase diagram for the quantum transitions as a function of tlie strength of the long-range 
interaction J'/Ji, for n = 4 and Ji = 1, J2 = 1, J3 = 1, J4 = 2. For J'/Ji > 0, the critical lines 
identify the first-order phase transitions and, for J'/Ji < 0, the second order phase transitions. 
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FIG. 7: Phase diagram for the quantum transitions as a function of tlie strength of the long-range 
interaction J'/Ji, for n = 4 and Ji = 1, J2 = 1, J3 = 2, J4 = 2. For J'/Ji > 0, the critical lines 
identify the first-order phase transitions and, for J'/Ji < 0, the second order phase transitions. 
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FIG. 8: (a) Magnetization as a function of h/Ji, and (b) universal curve for the magnetization 
as a function of the effective field hgff/Ji (heff = h + 2J'M^), at T = 0, for different values of 
r(r = J'/Ji), in the regions where the system undergoes first (r > 0) and second-order (r ^ 0) 
quantum transitions, for n = 4 and Ji = 1, J2 = 1, J3 = 1, J4 = 2. 
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FIG. 9: Phase diagram for the quantum transitions as a function of the strength of the long-range 
interaction J' / Ji, for n = 5 and Ji = 1, J2 = 1, J3 = 1,J4 = 1,J5 = 2. For J'/Ji > 0, the 
critical lines identify the first-order phase transitions and, for J'/Ji < 0, the second order phase 
transitions. 
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FIG. 10: (a) Magnetization as a function of /i/Ji, and (b) universal curve for the magnetization 
as a function of the effective field hgff/Ji {heff = h + 2J'M^), at T = 0, for different values of 
r{r = J'/Ji), in the regions where the system undergoes first (r > 0) and second-order (r ^ 0) 
quantum transitions, for n = 5 and Ji = 1, J2 = 1, J3 = 1, J4 = 1, J5 = 2. 
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h/J^ 

FIG. 11: Isothermal susceptibility Xt^^^ T = 0, as a function of h/Ji, for n = 3, Ji = 1, J2 = 1, 
J3 = 2, and different values of r(r = J'/Ji), (a) for r ^ and (b) for r ^ 0. 
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